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We discuss the recent Cabibbo's proposal to measure the tttt scattering lengths combination ao — «2 from the 
cusp effect in the 7r°7r° energy spectrum at threshold for K + — > 7r°7r°7r + and Kl — > 7r 7r 7r . We estimate the 
theoretical uncertainty of the ao — a2 determination at NLO in our approach and obtain that it is not smaller 
than 5 % for K + — > 7r°7r°7r + . One gets similar theoretical uncertainties if the neutral Kl — * 7r 7r 7r decay data 
below threshold are used instead. For this decay, there are very large theoretical uncertainties above threshold 
due to cancellations and data above threshold cannot be used to get the scattering lengths. 
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1. Introduction 

Final state interaction (FSI) phases at next-to- 
leading order (NLO) in Chiral Perturbation The- 
ory for K — > 37r are an important ingredient to 
obtain the charged CP-violating asymmetries at 
NLO [1-3]. The dominant contribution to these 
K + — > 3-7T FSI at NLO are from two-pion cuts 
and they were calculated analytically in [1]. 

Though to get the full K — > 37r amplitudes 
at order p 6 implies a two-loop calculation, one 
can get the FSI phases at NLO using the optical 
theorem within CHPT with the advantage that 
one just needs to know 7T7t scattering and if — > 3-7T 
both at 0(p 4 ). Notice that NLO in the dispersive 
part of the amplitude means one-loop and 0(p 4 ) 
in CHPT while NLO in the absorptive part of the 
amplitude means two-loops and 0(p 6 ) in CHPT. 

The study of FSI in K -> 3tt at NLO also be- 
came of relevance after the proposal by Cabibbo 
[4] to measure the combination ao — a 2 of 7T7t scat- 
tering lengths using the cusp effect in the 7r°7r° 
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spectrum at threshold in K + — > 7r°7r°7r + and 
Kl — > 7r 7r 7r decay rates. 3 Within this pro- 
posal, it has been recently presented in [6] the 
effects of FSI at NLO using formulas dictated by 
unitarity and analyticity approximated at second 
order in the 7T7t scattering lengths, et^ ~ 0.2. The 
error was therefore canonically assumed to be of 
order of af, i.e., 5%. There, they used a second 
order polynomial in the relevant final two-pion 
invariant energy S3 fitted to data to describe the 
K — > 37r vertex that enters in the formulas of the 
cusp effect. It is of interest to check this canon- 
ical error and provide a complementary analysis 
of this theoretical uncertainty estimate. 

In Ref. [7], we use our NLO in CHPT results 
for the real part of K — > 37r fitted to data to de- 
scribe the K — > 37r vertex that enters in formulas 
of the cusp effect. Notice that, as we want to ex- 
tract the pion scattering lengths, we do not want 
to predict the real part of K — > 3tt in CHPT at 
any order but to have the best possible descrip- 
tion fitted to data. We treat 7T7t scattering near 
threshold as in Cabibbo's original proposal. The 
advantage of using CHPT formulas for the fit to 
data of the real part of K — > 37r is that it con- 

3 The cusp effect in SU(2) tttt scattering was discussed in 
[5]. 
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tains the correct singularity structure at NLO in 
CHPT which can be systematically improved by 
going at higher orders. 

Contributions from next-to-next-to-leading or- 
der (NNLO) SU(3) CHPT in the isospin limit are 
expected typically to be around (3 ~ 5)%, so that 
our NLO results are just a first step in order to re- 
duce the theoretical error on the determination of 
the combination ao — a-i to the few per cent level. 
At NNLO, one can follow a procedure analogous 
to the one we use in [7] to get a more accurate 
measurement of ao — and check the assumed 
NNLO uncertainty. At that point, and in order 
to reach the few per cent level in the theoret- 
ical uncertainty, it will be necessary to include 
full isospin breaking effects at NLO too. These 
are also expected to be of a few per cent as was 
found for K 3tt in [8,9]. 

1.1. Basic notation for charged Kaon de- 
cays 

Near ir + ir~~ threshold, we can decompose the 
K + — > 7r°7r°7r + amplitude as follows [4,6] 



M)0+ 



( A 00+ + Boo+v±(s 3 ), for s 3 > 4m^ + 
1 A)o+ + iBoo + v±(s 3 ), for s 3 < 4ml + , 



(1) 



where Aqq + and i?oo+ arc in general singular 
functions except near the tt + tt~ threshold [10] 
and 
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Notice that these kinematical factors are taken 
with physical pion masses, in this way one can 
describe the cusp effect which is generated by the 
different behavior of K + — ► 7r°7r°7r + for the two 
neutral pions invariant energy above and below 
the S3 = 4to^ + threshold. 

With these definitions, the differential decay 
rate for this amplitude can be written as [6] 



|A) 0+ 1 2 = RcA 00+ + A A + v±(s 3 )A, 



CUSp > 



(3) 



with 



A 



cusp 



-2RcAoo+Im_Boo+ + 2ImA o+ReSoo+ , 
for S3 < 4ml+ ; 

2RcAoo+Re-Boo+ + 2Im J 4 o+ImBoo+ , 
for S3 > 4m?. + . 



A A = ImA 00+ +4(s 3 ) RcB 00+ + InxB 00+ 



The combination of real and imaginary ampli- 
tudes A cusp defined above parametrizes the cusp 
effect due to the ir + ir~ — > 7r°7r° re-scattering in 
the K + — ► 7r°7r°7r + decay rate. 

2. Discussion 

Our method is a variation of the original 
Cabibbo's proposal that uses NLO CHPT for 
the real part of K — > 37r vertex instead of the 
quadratic polynomial in S3 approximation used 
in [4,6] plus analyticity and unitarity. 

Notice that we do not use CHPT to predict 
the real part oi K — > 37r, but use its exact sin- 
gularity form at NLO in CHPT to fit it to data 
above threshold. If the two-loop CHPT singu- 
larity structure were known it could be used in 
order to take this singularity structure exactly in 
RcA o+- The treatment of nn scattering near 
threshold is independent of this choice and we 
treat it in the same way as in [6] . 

The cusp effect originates in the different con- 
tributions to K+ -> 7r 7r°7r+ and K L -> 7r 7r°7r° 
amplitudes above and below threshold of tt + tt~ 
production in the 7r°7r° pair invariant energy. We 
obtain these contributions using just analyticity 
and unitarity, in particular applying Cutkosky 
rules and the optical theorem above and below 
threshold to calculate the discontinuity across the 
physical cut. This allows us to separate irir scat- 
tering -which we want to measure- from the rest 
of K+ -» 7r°7r°7r+ or K L -> 7r°7r 7r°. 

The validity of the use of Cutkosky rules for 
in K — > 3ir decays is commented in [7]. In par- 
ticular, the real part of the discontinuity has a 
singularity when any of the Sj invariant energy 
reaches its pseudo-threshold at (mx — m 7r(i ) ) 2 as 
described in [10-12]. This singularity affects the 
description of the cusp effect using the disconti- 
nuity as dicussed in [7,10,11]. 

We would like to remark here that making the 
same approximations that were done in [6] we 
fully agree with their analytical results. In par- 
ticular, we checked that the use of the quadratic 
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polynomial in S3 in [6] produce negligible differ- 
ences -around 0.5 %- in A cusp in (4). 

We also pointed out that while the presence 
of that singularity at pseudo-thresholds does not 
affect Rc_Boo+(s3) and Re_B 000 (s3) when S3 is 
around threshold, one needs to take fully into ac- 
count its effects for <5ReAoo+ and <5ReA 000 when 
si or S2 is above {m 2 K — m\)/2. For a possible so- 
lution of this problem which does not simply use 
the discontinuity to describe the cusp effect see 
[10]. Another possibility could be to use, instead 
of 5Rcyloo+ and SRcA QQQ , the full- two loop (not 
available yet) finite relevant pieces to describe the 
singularities at thresholds at NLO in ReA o+ and 
RcA 000 , respectively. This could be fitted to data. 

In [7] we have also discussed the approxima- 
tions done in [6] and the numerical differences 
they induce in A cusp . We have found that though 
each one of them is individually negligible (be- 
tween 0.5 % to 1 %) they produce final differ- 
ences in the A cusp around 3 %. Of course, these 
approximations can be eliminated. 

Concerning the theoretical uncertainties in the 
determination of a — a 2 using our formulas, we 
concluded that for K + — > 7r°7r°7r + , this uncer- 
tainty is somewhat larger than 5 % if uncertain- 
ties are added quadratically and 7 % if added lin- 
early. I.e., we essentially agree with the estimate 
in [6] . Notice that we get our final theoretical un- 
certainty as the sum of several order 1% to 2% 
uncertainties to the canonical NNLO 5 % uncer- 
tainty. 

For the case Kl — > 7r 7r°7r°, we get that -if one 
uses just data below threshold- the uncertainty in 
the determination of a>o — a 2 is of the same order 
as for K + — ► 7r°7r°7r + . Above threshold, we found 
large numerical cancellations which preclude from 
using it. 

An expansion in the scattering lengths a» and 
Fcynman diagrams were used in [6] to do the 
power counting and obtain the cusp effect de- 
scription of K + — > 7r°7r°7r + at NLO. In general, 
when FSI tttt scattering effects are included at n- 
th order 4 , there appear new topologies in K — > 37r 
which give contributions to A cusp of order a™. 
The canonical uncertainty of the n-th order re- 

4 n = 1 order are for LO contributions. 



suits is thus af. Notice that the velocity factors 
that appear after applying the unitarity cuts can 
be order one -for instance v±((m K — m^)/2) ~ 
0.6 appear in Re i?oo+(4m^)- and do not suppress 
largely the naive af estimate. 

Our estimate for the uncertainty from NNLO, 
~ af, coincides numerically with the one made 
in [6], i.e. it is around 5 %. We conclude that 
one cannot expect to decrease this canonical 5% 
theoretical uncertainty of the NLO result unless 
one includes tttt scattering effects at NNLO. If 
one wants to reach the per cent level in the un- 
certainty of the determination of ao — a 2 from 
the cusp effect, one would need to include those 
NNLO re-scattering effects. As said above, at 
NNLO it is possible to follow a procedure analo- 
gous to the one we use here to get a more accurate 
measurement of ao — a 2 and check the estimated 
NNLO uncertainty. 

We have just included isospin breaking due to 
the different thresholds using two-pion physical 
phase spaces in the optical theorem and Cutkosky 
rules. This is needed to describe the cusp ef- 
fect. The rest of NLO isospin breaking is ex- 
pected to be important just at NNLO. At that 
order, isospin breaking effects in tttt scattering 
at threshold -both from quark masses and from 
clcctromagnetism- will have to be implemented 
and their uncertainties added. 

Finally, we believe that it is interesting to con- 
tinue investigating the proposal in [4,6] to mea- 
sure the non-perturbative tttt scattering lengths 
from the cusp effect in K + — > 7r°7r°7r + and Kl — ► 
7r 7r°7T . Another interesting direction is to de- 
velop an effective field theory in the scattering 
lengths which could both check the results in [6] 
and allow to go to NNLO. This type of studies 
is already underway and firsts results were pre- 
sented [10,11]. 
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